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Abstract 

Coccidiosis is a parasitic disease with major economic impact, one of whose main causative agents is Eimeria 
tenello. Chicl<en breeds display variable natural resistance to this disease. Unravelling the genetic bases of such 
variations could provide new clues for protection strategies. Transcriptomic experiments were conducted 
comparing resistant (Fayoumi) and susceptible (Leghorn) lines. Caecum and caecal tonsils were analysed. A global 
increase in differential gene expression following infection was observed for caecum comparisons, whereas a 
global decrease following infection was observed for caecal tonsils. 

Gene lists for infected tissues display 40 genes in common across breeds, 20 of which were specific to infected 
tissues. Among these specific genes, 9 belong to the 100 more differentially expressed genes of the infected 
caecum comparison. Gene expression networks were constructed in parallel, identifying highly connected genes. 
Comparing information from differential gene lists and gene network analysis allows one to highlight potential 
pivotal genes in the infection process, one of which was located in a putative significant QTL region for infection 
associated lesions. 



Background 

Coccidiosis is a disease with major global impact on sev- 
eral animal breeding species. It is one of the major dis- 
eases impacting poultry production with estimated of 
world losses averaging 1 billion $ / Year [1]. One of the 
main causative agents is the protozoan Eimeria tenella, 
an internal obligatory parasite, among other Eimeria 
species. During its life cycle, free stages invading epithe- 
lial cells alternate with a highly resistant stage (oocyst) 
adapted to long term environmental survival. The cae- 
cum is the primary site of infection. When infected, 
chickens {G alius gallus) present bloody droppings, stop 
feeding, face dehydration and huddle together. High 
mortality is mainly due to excessive blood losses, with 
recovering individuals likely to develop chronic illness. 
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Several therapeutic strategies are currently in use. 
Anticoccidial drugs exist but cases of drug resistance are 
developing, and their use is evermore prohibited by reg- 
ulations. Vaccination is partly possible, but it displays 
limitations due to multiple Eimeria species. Good farm- 
ing practices avoiding excessively dense animal concen- 
trations are highly recommended to this extent. 
Therefore, an improved knowledge of natural genetic 
resistance can pave the way towards protection against 
coccidiosis. Remarkably, chicken lines naturally display 
various degrees of resistance to E, tenella infection: the 
Egyptian Fayoumi line is known to be highly resistant 
compared to common layer breeds such as the White 
Leghorn [2]. Recent findings have located QTL regions 
related to lesions associated with this infection [3]. 
Here, we present a complementary approach using tran- 
scriptomic methods, which is the first to our knowledge 
to describe the effects of E. tenella infection on these 
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lines. These results strengthen and broaden previous 
results from one of the above mentioned QTL studies. 

Methods 

Experimental design 

Two tissues were studied: caecum and caecal tonsils (a 
lymphoid organ associated with caecum). For each tissue 
we first compared four to five Fayoumi (resistant line) 
and five White Leghorn (susceptible line) individuals 
that did not undergo infection, as control experiments. 
Then we performed similar comparisons for infected 
individuals at day 6 post infection. This design allowed 
us to make direct comparisons between lines, as well as 
between tissues and between infection status, and to dis- 
criminate between these different effects. 

Microarrays 

These transcriptomic experiments were performed using 
microarrays. Briefly, mRNA samples were directly 
extracted with Qiagen RNeasy Kit from disrupted tis- 
sues, including a DNAse step. mRNAs were checked for 
purity, quality and quantity prior to any further steps. 
Dye-swap comparisons were performed on chicken 20.6 
k oligo microarrays. Microarrays were analysed with an 
Agilent G2565BA Scanner at 10 |im resolution. Signal 
quantification was performed with the Genepix pro 6.0 
software. 

List of differentially expressed genes 

Data analysis was performed using the Bioconductor 
Limma package. Raw results were normalised as usual 
to overcome any cy5/cy3 fluorescence bias. Normaliza- 
tion was performed with the "normalizeWithinArrays" 
function with a loess correction for the dye effect and 
no background correction. No further normalization 
was required across arrays. The dye swap design implied 
technical replicates. Correlation structure between these 
technical replicates was accounted for using the "dupli- 
cateCorrelation" option. 

Gene networks 

Gene network reconstruction was based on Graphical 
Gaussian models and performed using the GeneNet R 
package [4]. 

Results 

Global profiling of differential expression 

In caecums, at a 5% Benjamini-Hochberg threshold, 
1473 genes were found to be differentially expressed 
between the two infected lines, and 174 genes were dif- 
ferentially expressed between the two non-infected con- 
trol lines. 42 of these genes were overlapping (present in 
both comparisons). 



In caecal tonsils, at a 5% Benjamini-Hochberg thresh- 
old, 61 genes were found to be differentially expressed 
between the two infected lines, and 113 genes were dif- 
ferentially expressed between the two non-infected con- 
trol lines. 17 of these genes were overlapping. 

Thus, the two tissues display a consistent remodelling 
of their differential gene expression during infection, as 
opposed to the basal differential gene expression which 
is of a similar magnitude (174 and 113 genes) in the 
two tissues. 

However, the two responses display apparently oppo- 
site trends: a global increase of differentially expressed 
genes is observed (from 174 to 1473 genes) in caecums, 
whereas a global decrease is observed (from 113 to 61 
genes) in caecal tonsils. This might be related to the fact 
that the two tissues are not impacted at the same time 
of the infection process, and will be further investigated. 
But due to the quite low overlap of genes between 
infected and control experiments, the main and com- 
mon feature of these infection effects remains the global 
remodelling of differentially expressed gene list. 

In order to narrow the core group of genes globally 
involved in the infection process, we compared the dif- 
ferent gene lists. 40 genes were found common to both 
tissues' "infected" comparisons. 20 of these genes were 
only found in both "infected" comparison lists (that is to 
say not in the "control" lists), and therefore can be 
labelled as "specific". Among these specific genes, 9 
belonged to the 100 most differentially expressed genes 
of the "infected" caecum comparison, and therefore can 
be labelled as "top specific". 

First insights in gene network analysis 

Networks were obtained using the 100 most differen- 
tially expressed genes between the two infected lines for 
the caecum. Significant edges were found at a 10% and 
20% local false discovery rate (FDR) threshold. The cor- 
responding networks (displaying respectively 8 and 18 
genes) are presented in figure 1. 

Comparing gene interaction network and differential 
expression information 

4 of the above described 9 top specific genes were found 
in the 20% network, and their corresponding proteins 
are: HPSE (GGA4 - heparan diving protein involved in 
leucocytes migration pathway), SORF2 (GGA18 - maps 
in a highly significant QTL region for lesion [3]), 
ABCB4 (GGA2 - defence against xenobiotics), DUSP26 
(GG A22 - signalling pathway modulation). In addition, 2 
of the 40 common genes were found in the 10% net- 
work: TECTB (GGA6 - matrix protein) and an ECR 
(GGA2- evolutionary conserved region/no gene annota- 
tion- similar to a replication licensing factor). 
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Figure 1 Gene Expression networlcs of the caecum infection by E tenella. (A) Significant edges at 20% FDR (B) Significant edges at 10% 
FDR. Plain lines: positive relationship; dashed lines: negative relationship 



Conclusions 

These first results show the usefulness of performing 
and integrating different transcriptomic methods when 
studying the genetic bases of natural resistance to dis- 
eases. Microarrays allowed a global view of transcrip- 
tomic modification upon infection, detecting different 
trends and amplitude of gene expression changes in the 
two tissues. As a starting point, we focused on the com- 
mon genes differentially expressed in both tissues. A 
detailed analysis of the differentially expressed genes 
lists in each tissue, including a clustering based on varia- 
tion profiles, is currently ongoing. Comparing this infor- 
mation, enhanced by a two-tissue comparison, with 
gene network analysis obtained with the same list 
allowed the identification of core genes significantly 
associated with the day 6 post infection molecular phe- 
notype. It is noteworthy that several of these genes are 
related to extracellular matrix (hereby not excluding leu- 
cocyte migration pathways associated with infection), 
and one of them, SORF2, maps into a previously 
described highly significant QTL region for lesions. This 
can be seen as a partial validation of these previous 
results, and demonstrates that combined QTL/transcrip- 
tome studies can lead to identification of candidate 
genes. Nonetheless, further refinements of these results 
are required, in order to better understand the biological 
meaning of the obtained networks. The status of highly 
differentially expressed genes not found in the gene net- 
works will need special further investigation. In parallel. 



an extensive comparison with the above mentioned 
QTL data will help to understand why some regions are 
congruent in the two types of studies, and others not. In 
depth analysis of the gene network relationships will 
bring a more subtle view of the genome dynamics corre- 
lated with resistance or susceptibility to coccidiosis. 
Transcriptomics of intermediary steps of the infection 
kinetics will refine these initial results, provided that 
mRNA expression can be observed preceding clinical 
signs of infection such as lesions. Our results describe 
infection by Eimeria tenella, however other Eimeria spe- 
cies are alternative causative agents of coccidiosis. Com- 
paring and interpreting similar results obtained with 
these different species might help decipher the specific 
and common features of coccidiosis infection at the 
molecular level. 

These first results demonstrate the efficiency of a 
microarray strategy to identify genes and networks 
involved in infection. Further, they underline the useful- 
ness of local breeds as genetic resources. 
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